Preferred States of Decoherence under Intermediate System-Environment Coupling 
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The notion that decoherence rapidly reduces a superposition state to an incoherent mixture im- 
plicitly adopts a special representation, namely, the representation of preferred (pointer) states (PS). 
For weak or strong system-environment coupling, the behavior of PS is well known. Via a simple 
dynamical model that simulates a two-level system interacting with few other degrees of freedom as 
its environment, it is shown that even for intermediate system-environment coupling, approximate 
PS may still emerge from the coherent quantum dynamics of the whole system in the absence of 
any thermal averaging. The found PS can also continuously deform to expected limits for weak or 
strong system-environment coupling. Computational results are also qualitatively explained. The 
findings should be useful towards further understanding of decoherence and quantum thermalization 
processes. 
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Introduction — As illustrated by the Schrodinger cat 
paradox, there is a clash between the quantum superpo- 
sition principle and the way we perceive the macroscopic 
reality. So how does a classical world emerge from the 
quantum? One promising solution to this profound ques- 
tion is decoherence, i.e., the loss of quantum coherence 
due to the interaction of a system of interest with its 
environment [IH3- Decoherence may rapidly reduce a 
coherent superposition state of the system to an incoher- 
ent mixture. During this process the environment singles 
out special basis states, often called "preferred (pointer) 
states" (PS), of which a classical probabilistic description 
becomes sufficient to describe the system and the bizarre 
superposition state of the PS is out of the picture. That 
is, in the PS representation the reduced density matrix 
(RDM) of the system becomes diagonal as time evolves. 

Such a decoherence perspective is not expected to re- 
solve all conceptual issues regarding quantum weirdness 
vs classical reality. Nevertheless, it is highly useful as it 
implies the environmental-dependence of the quantum- 
classical transition and the representation-dependent na- 
ture of decoherence issues. For example, different envi- 
ronments may select different PS, and a quantum system 
decohered in one PS representation may still possess cer- 
tain quantum coherence in other representations. Going 
further, one may envision the possibilities of environment 
engineering to form desired PS such that system prop- 
erties are robust to decoherence. 

PS have been identified in several cases. If the system- 
environment interaction is in the adiabatic limit or if it 
commutes with the system's self-Hamiltonian, then pop- 



ulations on the energy eigenstates of the system's self- 
Hamiltonian do not change but their relative phases are 
destroyed by the environment. The energy eigenstates 
then form the PS [oT-Hoj]. Analogous to this, if the system- 
environment coupling is sufficiently weak, then the en- 
ergy relaxation time scale is typically much longer than 
the pure-dephasing time scale. As a result, the energy 
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eigenstates still form the PS before relaxation sets in [11 
1 1 31 ] . Again, for weak system-environment interaction but 
for a longer time scale, a model of quantum Brownian 
motion reveals that coherent states localized in both po- 
sition and momentum turn out to be the PS [H,[2, 111 14 1- 
In the opposite situation, the system-environment cou- 
pling is strong and the system's self-Hamiltonian becomes 
negligible within a certain time period. In this case, the 
eigenstates of the system-environment interaction Hamil- 
tonian form the approximate PS 0-0] • 

Little has been said about the possible existence of 
PS for a generic system-environment coupling (i.e., not 
commutable with system's self-Hamiltonian) of interme- 
diate strength. Under such a situation, the widely used 
quantum-master-equation approach or other perturba- 
tive approaches may not be applicable in analyzing the 
existence of PS. Also motivated by the ongoing investi- 
gations of quantum thermalization processes 15|, [16[ , we 
choose to work with a simple dynamical model to ad- 
dress the issue of PS. That is, within a single isolated 
quantum system composed of interacting quantum sub- 
systems, will the concept of PS still work well in describ- 
ing the decoherence of one subsystem due to its interac- 
tion with other subsystems 1 7jj ? 

We start from a computationally intuitive definition 
of PS. We then show interesting evidence that PS may 
still exist for intermediate system-environment coupling 
and further explain why this is possible. The found PS, 
neither the system's energy eigenstates nor the eigen- 
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states of the system-environment interaction Hamilto- 
nian, undergo continuous deformation as the system- 
environment coupling strength varies. These findings 
show that decoherence-induced superselection rule can 
be twofold: superposition states of the PS are destroyed 
but PS themselves can be rich superposition states for 
intermediate system-environment coupling. Equally in- 
teresting, it can be concluded that the concept of PS 
is still important in understanding quantum dynamical 
processes in the absence of any thermal averaging. 

Identifying PS from the time-evolving RDM — If PS 
exists, then the RDM will gradually become diagonal in 
the PS representation. On the other hand, the same 
RDM is always diagonal in its own eigenrepresentation. 
Therefore, if we computationally track the eigenstates 
of the RDM, then we can see clearly whether or not a 
well-defined set of PS can emerge from a decoherence 
process. That is, if after a certain period the eigenstates 
of the RDM are found to evolve closely around a fixed 
basis set, then this fixed set of states can be defined as 
the PS, at least approximately. This computational def- 
inition of PS extends the PS criterion used by Diosi and 
Kiefer [8]. Note also that such a definition of PS gives 
up their precise analytical form. Consistent with this 
picture, the off-diagonal elements of the RDM in the PS 
representation must be also small when compared with 
the difference of its diagonal elements. Indeed, were the 
RDM diagonal elements almost degenerate, then a small 
fluctuation in the off-diagonal elements can still cause 
a drastic rotation of the RDM eigenstates, a fact that 
would contradict with the existence of PS 18| ■ We hence 
mainly work in the parameter regimes where an appre- 
ciable difference between the diagonal elements of RDM 
can emerge from the dynamics. These preliminaries also 
make it clear that even a stationary RDM does not nec- 
essarily mean the existence of PS. 

Consider now a two-level system S interacting with its 
environment £, with a total Hamiltonian H = Hs + Hi + 
Hg, where H$ and Hs are the Hamiltonians of S and £, 
and Hi is the system-environment interaction Hamilto- 
nian, with [Hs,Hi] ^ 0. Eigenstates of H$ are denoted 
by |-Efc), with k = 0, 1. Throughout, we use \$>) to denote 
a state vector for the whole system-environment combi- 
nation, denoted by S + £ and isolated from any ther- 
mal bath. The time-evolving RDM for S is given by 
p s (t) = Tr £ (|*(t))(*(t)|), where = e - im / H \^(0)) 

coherently evolves according to Schrodinger equation for 
S + £. 

Eigenstates of the RDM p s (t) are represented by 
\pk(t)}, with eigenvalues pk(t), i.e., p s {t)\pk(t)) = 
Pk(t)\pk{t)}, with k = 0, 1. The two states \po(t)) and 
\pi(t)) also form an orthonormal basis set for the Hilbert 
subspace associated with Hs- The distance between this 
basis set and another basis set \r)k') (k' = 0,1) for the 
same subspace may be measured by D(\pk{t)),\r]k')) = 
1 — \(pk(t) \Vk'}\ i with k and k' determined by the con- 



dition \{pk\Vk')\ 2 > 1/2. A time-averaged distance d 
over a period [£ a ,tf,] can then be defined as d(\rjk')) = 
[1/(4 - ta)]j£j& D(\p k (t')),\r) k .}). If, for a particu- 
lar basis set {|/5/c')} defined below, d(\pk')) is small for 
sufficiently large t a , and % — t a , then p s (t) becomes 
almost diagonal in the {|po), \pi)} representation, and 
hence {|/5o)> I Pi}} can be computationally identified as 
the PS. 

To find {\po), Pi}} that may reflect the average behav- 
ior of \pk(t)) with acceptable fluctuations, we calculate 
the time-evolving RDM eigenstates \pk{t)), average the 
density matrix \pk(t)){pk{t)\ over time (value of k = or 
k = 1 is chosen to maintain a continuity), and then obtain 
a time-averaged density matrix p. Finally, the eigenstates 
of p are defined as the basis states {\po),\pi}} [19J- If, in 
the {\p ), |pi)} representation, \p s 01 (t)\ < |p5o(*) ~ Pii (*) I 
for sufficiently large times t, then d(\pk')) is small and PS 
can hence be identified [3]. If this is not the case, then 
PS fails to emerge from the dynamics. 

Model — We now turn to a concrete model. To reflect 
the fact that typically a small system S is not directly 
coupled to the whole of its environment £ , we let S be 
directly coupled to a small component A of £, and then 
let A be further coupled to the rest part B of £ , with 
£ = A + B. For convenience, A is also assumed to be 
a two-level system. Such kind of coupling scheme was 
also considered recently [2(| to model a nonlinear system- 
environment coupling. It can be also qualitatively argued 
that our coupling scheme can yield much less fluctuation 
in the RDM than a full coupling between S and £ does. 
The B part of £ is simulated by a quantum kicked rotor 
on a torus with only 1 degree of freedom, whose classical 
" 22 



limit is fully chaotic [21 



The irregular motion of B 
due to quantum chaos, instead of many noninteracting 
degrees of freedom of a thermal bath, is responsible for 
decoherence in S. In terms of standard Pauli matrices 
and operators for a kicked rotor in dimensionless units, 
the Hamiltonians for the system, the environment, and 
their coupling are 



H s = uj x a* +u z a b zl Hi = ea 
H e = H A +H B +H AB . 



(1) 



Here, Ha = w^af, H B = \ + v cos 7 ^ 6(t - jT), and 
Hab — cos 7 J] j 8(t — jT), where p and 7 are mo- 
mentum and coordinate operators of the kicked rotor. 
Since the system-environment coupling is already of the 
erf type, for generality Hs is made to contain both erf 
and erf terms [23]. The unitary propagator associated 
with one period T is (with h = 1). 



jj t _ e -iT(w a fff+u(,trf +WA<r^+ecrf &<r£) 
iT^r „—iv cos 7 — cos 7 



xe 



2 e 



(2) 



The initial state is chosen as 1^(0)) = \iPq) ( E)\0)a < S)\ ( Po)b, 
where |-0f ) and |<po)s are vectors in the Hilbert spaces of 
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FIG. 1: (Color online) Angle # (in unit of radiant) to be ro- 
tated to reach state |pi) (empty symbols) from one eigenstate 
of Hs (square) or of Hi (triangle), for a wide range of e. 
State |pi) is one numerically found eigenstate of p (a time- 
averaged RDM for t € [30000T, 40000T]). For comparison, 
the 9 values to reach the state \a) (filled symbols) determined 
theoretically [details after Eq. (6)] by maximizing ||Aif||, are 
also presented. The initial state of the system S is placed 
in a superposition state 0.8exp(5i)|s+) +0.6|x_), where \x+) 
and \xJ) represent spin-up and spin-down states along the 
a;-direction. For the initial state of the environment, A is 
placed in an eigenstate of Ha and the kicked-rotor state is cho- 
sen randomly. Other system parameters are uj x = 0.5 x 10 3 , 
oj z = 1.0 x 10 3 , cja = 1.5 x 10 3 , v = 90/T,N = 2 12 , A = 0.1. 
Inset: Decay of the off-diagonal element of RDM with time, 
in \Ek) representation (upper red curve), in eigenrepresenta- 
tion of Hi (middle blue curve), and in representation of the 
PS identified here (bottom dotted curve), for e = 2000. Note 
that the decay is essentially done within about 600T. 



S and of B, and \0)a is an eigenstate of Ha- The quan- 
tum kicked rotor is quantized on a phase space torus 
with a Hilbert space dimension N = 4096, whose initial 
state is taken as a randomly generated vector from its 
Hilbert space (quantum recurrence time is already suffi- 
ciently large). The kicking period T is taken as 2ir/N. 



Typical values of u> x 



and uia are set around 10 in 



dimensionless units, such that within one kicking period 
T the characteristic phase evolution of the two-level sys- 
tems are of the order of unity. Many initial states were 
studied but in Fig. 1 we only report representative results 
for one initial state. Note also that one key parameter is 
e, which represents the strength of system-environment 
coupling. 

As shown in Rcf. 



if the system-environment cou- 



pling strength e is below a threshold e p , then the off- 
diagonal elements of the RDM in the eigenrepresentation 
of Hs will show a Gaussian-type decay. The dephas- 
ing time of energy eigenstates (T 2 ) then scales as e" 1 , 
whereas the population relaxation time (Ti) goes as e~ 2 
(obtained from Fermi's golden rule). In our model we 
find e p ~ l/N. So for e < 1/N, T 2 <C T u and hence 
the energy eigenstates \Ek) form the PS. Detailed cal- 



FIG. 2: (Color online) Upper panel: Behavior of the time- 
evolving eigenstates of the system's RDM, as described by 
a time-averaged distance d from the states (|po), |pi)) found 
computationally in Fig. 1 (circles), from eigenstates of Hs 
(dashed red line), and from eigenstates of Hi (dotted blue 
line), for a wide range of e. Other system parameters are 
the same as in Fig.l. Bottom panel: distance d between 
time-evolving eigenstates of the system's RDM and states 
(|po), |pi ) ) found computationally in the case of A = 1.0. 
Large d values in the bottom panel indicate the loss of PS. 



culations from our present model confirm this and also 
reveal something interesting. As shown in Fig. 1 (empty 
squares), \E^) are found to agree well with the compu- 
tationally found PS (i.e., very small values of 6*) for e 
as large as 10 1 . Though our previous work [ll[ did not 
rule out the possibility of \Ek) being the PS for e > e p , 
it is remarkable to see that \Ek) here still form the PS 
even for e ^> e p . This should be related to the fact that 
here the system (S) is only directly coupled with a small 
component (A) of the environment. 

Results for larger values of £ are also detailed in Fig. 1. 
Consider first e values approaching 10 4 , i.e., the right end 
of the curve shown with empty triangles. In these cases, 
the computationally found states \pk) are rotated from 
the Hi eigenstates by essentially a zero angle. Hence, 
the eigenstates of Hi can be regarded as the PS in this 
strong coupling case, even for a time scale much larger 
than the characteristic scale of Hs- Next we turn to 
intermediate cases with e E [10 2 ,5 x 10 3 ]. As seen from 
Fig. 1, states \pk) can notably deviate from \Ek) as well 
as the Hi eigenstates. As we tune up the value of e, 
states \pk) exhibit a clear and smooth transition from 
being close to \Ek) to being close to the Hi eigenstates. 
Further, as a consistency check, the inset of Fig. 1 shows 
the decay of the off-diagonal elements of the RDM in 
three representations, for e = 2000 as an example. It is 
seen that only in the \pk) representation, the off-diagonal 
elements decay to small values with some fluctuations 
24|. 
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It is yet to be shown that at sufficiently later times 
the eigenstates \pk(t)) of the time-evolving RDM only 
slightly fluctuate around |p/-). The upper panel of Fig. 
2 depicts the distance d(\p~k>)) vs e (solid line) (i.e., for 
the PS identified in Fig. 1). It is seen that for the en- 
tire considered regime of e, d remains impressively small. 
For intermediate values of e, it is much smaller than the 
same d-distancc between \pk(t)) and eigenstates of Hs 
(dashed red line) or between |p/s(t)) and eigenstates of 
Hj (dotted blue line). Therefore, after an initial period 
the time-evolving RDM eigenstates {\po(t)),\pi{t))} do 
remain close to {\po), \pi)}, suggesting that the RDM be- 
comes almost diagonal in the \pk) representation, even for 
intermediate values of e. We finally infer that the states 
\pk), whose behavior is shown in Fig.l, are indeed excel- 
lent PS emerging from the coherent quantum dynamics 
of S + £ . Other detailed calculations also indicate that 
PS does not always exist. For example, if the coupling 
between A and B is also very strong (e.g., A = 1.0, bot- 
tom panel of Fig. 2), then the diagonal elements of RDM 
become very close, d can reach as large as 0.25 (dashed 
line) (a value that can be estimated theoretically 18]), 
and consequently PS is lost. 

Theoretical insights — We shall now develop some in- 
sights into our computational results. We rewrite the 
total state for S + £ as 

\V(t)} = \a}\<p a (t)) + \f3)\Mt)), (3) 

where (\a),\f3)) is a chosen time- independent orthonor- 
mal basis set in the Hilbert subspace for S, \<fi a (t)) and 
\4>p it)) are the associated "expansion states" living the 
Hilbert subspace for £ . The time-dependence of the off- 
diagonal element of the system's RDM is then given by 

p s a p = (a\p s \P) = (Mt)\<i>«(t)). (4) 

With this notation, seeking PS then becomes the search 
for (| a), |/3)), such that the evolution of \4> a {t)) is as dif- 
ferent as possible from \<f>p(t)}. This insight motivates 
us to examine the time evolution of the states \(f> a (t)) 
and \<f>p{t)) used in Eq. ©. To that end, we first define 
HS p = (a\H s \P), Hip ee (a\H!\p) (and those by a O /3). 
Note that H^p thus defined is a scalar, but H^p is still 
an operator on the Hilbert subspace for £. Just to have 
a rather compact Schrbdinger-like equation for \(f) a (t)) 
and \4>p(t)), we introduce more operators on the £ sub- 
space, i.e., H aa = H £ + Hl a + H% a , H a p = H^ + H^p, 
K a ee H a pHppH~p, and J a = H a pHp a . Using these 
definitions and the Schrodinger equation for S + £ , we 
obtain 

*4-|<M = H aa \4> a ) + i\£a), (5) 
dt 

where |£ a ) = —iH a p\(f)p), with 

i-lllta) = Kalta) -iJa\<t>a)- (6) 



An analogous equation for \<fip(t)) is obtained by exchang- 
ing a and /3. 

Equation ([7]) indicates that the difference between the 
evolution of \(/> a {t)) and that of \<j)p{t)) is caused by the 
difference between the operators H aa and Hpp and by 
the difference between |£ a ) and |^). Further, Eq. (j8]) 
shows that |£ a ) and \£p) evolve differently due to the 
difference between K a and Kp and between J a and 
J p. To quantify these operator differences we define 
AH ee H aa - Hpp, AK = K a - Kp, and A J = J a - Jp. 
In our model, because Hj is a direct product of two spin 
operators, one finds A J = Furthermore, for s 

much less or much larger than the energy scale of Hs, 
we find AK w —AH [18|. For cases with intermediate 
values of e, a more detailed analysis (3) gives again that 
AK w —AH, at least for those states (la), \(3)) that max- 
imize ||Ai?||, where ||A/f|| represents the Frobenius-2 
norm, a simple measure of AH. Putting these observa- 
tions together, we intuitively expect (not a proof) that 
for the entire considered regime of e, the basis states 
(|a), |/3}) maximizing ||Ai?|| may approximately give the 
most substantial difference between \<j> a {t)) and \<j>p{t)) 
and hence the most significant decay of \p s a p\- As such, 
the basis states (|a), |/3)) theoretically determined by 
maximizing ||AiJ|| should agree with the PS computa- 
tionally obtained above. 

In Fig.l we compare the PS (empty symbols) found 
from the decoherence dynamics with the states (|a), |/3)) 
(filled symbols) determined directly by maximizing 
||AJJ||. In terms of their relation with the eigenstates of 
Hs and of Hj, nice agreement is obtained for the whole 
regime of e. 

It is also interesting to discuss the implication of the 
term i\£ a ) m Eq. (|7J|. Note that this term does not 
preserve the norm (4> a \4> a )- Hence, populations on the 
basis states (|a), |/3}), even when they are identified as 
the PS, can still fluctuate with time. This constitutes a 
crucial difference from a pure-dephasing problem. Un- 
like in a pure-dephasing problem, here the decay of 
Pap = (^Wl^aW) cannot be interpreted as that of the 
overlap of two independent environment histories \<f> a (t)) 
and \<f>p(t)). Instead, these two evolution histories are 
mingled together through population transitions between 
them. Difference from a pure-dephasing picture is also 
made evident by the role of the term AH S = H^ a — Hpp 
as one component of AH = AH S + H^ a — Hip. For a 
pure-dephasing problem, i.e., if the term «|£ a ) is switched 
off, then the component AH S becomes irrelevant: it 
is a c-number for the environmental Hilbert subspace, 
and hence cannot cause the decay of \{4>pit)\(j) a it))\. By 
contrast, in our model the term AH S in AH is found 
to be necessary for predicting the PS with intermedi- 
ate system-environment coupling. That is, without this 
term, maximization of ||Ai?|| would incorrectly predict 
that eigenstates of Hi are the PS regardless of the value 
of e. 
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Conclusion — The concept of PS may still apply if an 
environment with many degrees of freedom (like a ther- 
mal bath) is replaced by an environment with very few 
degrees of freedom. Approximate PS are shown to exist 
for intermediate system-environment coupling and can 
continuously deform to expected limits. Such types of 
PS emerging from quantum dynamics alone are of impor- 
tance to understanding decoherence and thermalization 
processes. 
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Supplementary Materials I: Some operators and 
equations used in the manuscript 

To qualitatively understand the emergence of preferred 
(pointer) states (PS) for intermediate strength of system- 
environment coupling, we first recapitulate some equa- 
tions and operators defined in our manuscript. In partic- 
ular, the "expansion states" living in the Hilbert sub- 
space of the environment is denoted by \4> a (t)), with 
|¥(f)) = \a)\<f> a (t)) + \P)\Mt)), where (\a), |/3)) is an 
time-independent orthonormal basis set in the Hilbert 
subspace of the system S. 

Using Schrodinger equation for the system and the en- 
vironment as a whole, under the condition that H a p does 
not explicitly depend on the time t and is a reversible op- 
erator, the equation of motion of \<p a {t)) is given by 

i^\4>a) = H aa \4> a ) + i\£ a ), (7) 

i^\U) = K a \U)-iJ a \<t> a ), (8) 

where 

H aa = H £ + H^a + H aa , 
K a = H a pHppH~p, 

Ja = H a pHf3 a , (9) 

with 

H a f3 = H a0 + H a p, 
Hip = (a\H s \p), 

Hi p = (a\Hi\p), (10) 

and analogous expressions by exchanging a and /3. In 
the manuscript we then defined the following operators: 

AH = H aa — Hp , 
AK = K a — Kp, 

AJ = J a - J . (11) 
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The vector \£ a ) is in fact related to \<j>p) by |£ Q ) = 
—iH a p\<j)p). For an initial system-environment direct- 
product state |*(0)) = |^o)l^o)) w e have |0 O (O)) = 
{(a\^))\4>o) and |^(0)) = ((/3|^))l0o). Thus, the ini- 
tial condition of Eqs. {7} and ([SJ is given by |0 a (O)) and 

|ea(0)) = -i((/3|^»^l<^0). 

Note that (<f> a (t)\(j) a (t)) and give the 

probabilities of finding the system in the states \a) and 
|/3), respectively. During the time evolution, these two 
quantities need not to be conserved because the problem 
of interest here is not a pure-dephasing problem. 



Supplementary Materials II: On conditions for PS to 
be identified computationally 

In our manuscript we also mentioned conditions under 
which we can computationally identify PS from the time- 
evolving reduced-density matrix (RDM) of the system. 
These conditions are further elaborated in this section. 
After some simple derivation, we obtain the following ex- 
pression for the distance D(\pk(t)), \ a)) = 1 — \(pk(t)\a)\ 2 , 
where the eigenstate \pk(t)) of the RDM p s (t) is chosen 
by the condition \{pk(t)\a)\ 2 > 1/2, 



D(\p k (t)),\a)) = 



5p+^{5p) 2 +i\p a aP \ 2 \ +\p s aP \ 



.(12) 



Here, Sp = \p* aa - p s P p\, with p s aa (p s a0 ) being the diagonal 
(off-diagonal) matrix elements of the RDM. The time- 
averaged value of D(\pk(t)}, \a)) is denoted by d(\a)) in 
the manuscript. 

Suppose at long times, the overlap ((f>p(t)\<f> a (t)) is al- 
ready of a small magnitude. Then, p s a/3 = (cf)p(t)\4> a (t)) 
is small by definition. Furthermore, we assume that the 
populations in the two basis states, namely, p s aa and p^, 
have also approached some equilibrium values at long 
times and can maintain an appreciable difference Sp such 
that Sp > \p s aa \. Then, Eq. (IS]) becomes 



D(\ Pk (t)),\a)) 



\Pap[ 

(s P y 



This then gives 



d(\a)) 



1 



tb — t a 



dt 



(s P y 



< i. 



(13) 



(14) 



Sufficiently small value of d(\a)) indicates that eigen- 
states of the RDM p s (t) only have small fluctuations 
around the fixed set of basis (|a), |/3}). This fixed set 
of basis states then form the PS, at least approximately. 
In our manuscript we computationally identified such a 
basis set, which is denoted by \pk)- 

It is now clear that PS can be computationally identi- 
fied under the following conditions: (1) small \p s a a\ and 



(2) relatively large Sp. Smallness of \p s a o\ is equivalent to 
the smallness of the overlap \((f)p(t)\(f) a (t))\. This occurs 
if AH, AK, and AJ are made as different as possible, 
because these operators will determine the evolution dif- 
ference between \<j)p(t)) and \4> a (t)). 



Supplementary Materials IILDistance d averaged 
over the Hilbert space 

In the case of the bottom panel of Fig. 2, the envi- 
ronment part A is also strongly coupled with the kicked 
rotor, the distance d is large and hence the PS is lost. 
Here we give a simple explanation. Due to the strong 
coupling within the environment and the strong cou- 
pling between S and £ , we can assume that the mo- 
tion is sufficiently random in the total Hilbert space for 
S + £ . Therefore, the long-time-averaged d can be ap- 
proximated by a Hilbert space average. Let (d) denote 
this total-Hilbert-space-averaged distance describing the 
difference between RDM eigentates and a fixed set of ba- 
sis states (\a), \/3)). In the following, we estimate (d) us- 
ing a method similar to the canonical typicality approach 
discussed in Ref . [l| . The dimension of the Hilbert sub- 
space of the environment £ is denoted by ng and the 
dimension of S + £ is 2ri£ . 

Assuming that the dynamics is almost ergodic on the 
total Hilbert space, a "typical" vector in the total Hilbert 
space can be written as 

\*)=tf y £\a)\$i), with|^)=^C aj -|i), (15) 

a j 

where {\j}} is an orthonormal basis set in the Hilbert 
subspace for the environment £ , j\f is the normalization 
coefficient, 

(16) 



j\f - 



J2aj \Caj\ 2 



2n £ 



and the real and imaginary parts of C a j are assumed 
to be independent real Gaussian random variables with 
mean zero and variance 0.5. For this typical vector, the 
elements of the RDM are given by 



PaP 



Sp 



\Paa Ppp\ 



■ jV 2 



\a 



Ay I 



(17) 



(18) 



In the limit of large ng , after simple derivation one finds 

Ip^I 2 ) 



lim 

ne — >oc 



(Sp 2 ) 



1 

2' 



(19) 



with the average (•) taken over all typical vectors defined 
above. Furthermore, using Eq. (|12l) . one obtains that 
(d) = 0.25, a result consistent with the bottom panel of 
Fig. 2. This clearly shows that PS does not always exist. 
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Supplementary Materials IV: Analysis of AH, AK, 
and A J 



Finally, we discuss how the three difference operators 
AH, AK, and AJ may be used to predict PS without 
actually following the decoherence dynamics computa- 
tionally. The qualitative picture is that in the PS rep- 
resentation, the operators AH, AK, and AJ should be 
maximized such that the two evolution histories \4>p{t)) 
and \4>a(t)) are as different as possible and hence the off- 
diagonal elements of the RDM will decay to very small 
values. Because in our model, Hj has a product form, 
namely, Hj — eerf ® a^, we have 

[H^Hip] = 0, (20) 
A J = H aSi Hj ia - HjiaHaji = 0. (21) 

Given A J = 0, we focus only on the maximization of AH 
and AK to understand the emergence of PS. 

Interestingly, the two operators AH and AK are also 
related to each other. Let X a — \H a p,Hpp], we find (af- 
ter some straightforward calculations) the following ex- 
pression 



where 



AA = -AH + Y, 



Y - X a H a l - X m 



(22) 



(23) 



For our model, we may further simplify X a and obtain 
X a = [H aj3 , Hg}. Next, to have explicit expressions of op- 
erators, we now take the eigenstates of a x as our working 
representation, where \ l) x and \0) x represent spin up and 
down along the x-axis. We expand a given set of basis 
states { | a), as the following: 



|a> = ae^|l) x + 6|0) x , 
|/3) = &e^|l) x -a|0) x , 



(24) 



where a, b G [0, 1] are real expansion coefficients. Us- 
ing these two explicit expressions and our direct product 
form of Hi, we find 



X a = [(l- 2b 2 ) sin(p) + i cos(</?)] {2ieu} A )<Jy (25) 



and 



Y = 2isu A (i]L — f]L), 



(26) 



where <jy is the Pauli operator along the y-axis for the 
two- level system A as a part of the environment, 

ci _ c 2 



V = 



d\ 



£ 2 C? ' V d% 



2„2 ■ 



ci = (1 — 2b 2 ) sin(t^) + i cos(<p), 
di = 2uj x \/l - b 2 b + u z ci, 
c 2 = c\, d 2 = d\, 



(27) 




FIG. 3: Variation of the ratio r with e for an arbitrarily 
chosen basis (|a),|/3)). Parameters: the same as in Fig.l of 
the text. 



and 



-iecia x + d\(jy, L ■■ 



-iec 2 a + d 2 a 



v ' 



(28) 



Let us consider the first limiting case of e <C u x . Ne- 
glecting terms of the order of e/uj x and noticing that u>j 
(i = x, z, A) are of the same order of magnitude, we find 



n 



Cl 



£2 
d 2 



L ~ d\o , L ~ d 2 a J 



(29) 



Thus, 



Y ~ 2ie 



( uj a ci 
V di 



U A C 2 \ A 



e. 



(30) 

uj x , one infers that A A' ~ 



Clearly then, because AH 
-AH. 

In the second limiting case of e ^> uj x , we neglect all 
terms of the order of oo x /e and then obtain 



V 



1 



9 

e z ci 



V ■ 



e 2 c 2 



-iscia x 



L 



A 

-i£c 2 a x . 



Thus, 



Y ~ 2lo a o^ - 2uj a ct^ = 0. 



(31) 



(32) 



This also leads us to the expectation that AA ~ —AH. 
That is, for very small or very large system-environment 
coupling strength e, we always have AA ~ —AH. 

The situation of intermediate system-environment cou- 
pling strength with e ~ co x is more complicated. To bet- 
ter understand this case we have performed numerical 
calculations of the ratio, 



r(e,|a» 



11*1 
\AH\ 



(33) 
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FIG. 4: Variation of the ratio r with e for the basis (|a), \/3)) 
that maximizes the norm ||A_ff||. 

where || . . . | indicates a Frobenius-2 norm. Figure 1 
shows the dependence of the ratio r on e for an arbi- 
trarily chosen basis (|a), It shows that the ratio can 
be large for intermediate values of e. Therefore, AK may 
have considerable deviations from —AH for a general ba- 
sis set (\a), |/3)). On the other hand, for the basis set that 
maximizes the norm ||Aff||, it is found that the ratio r 
is negligibly small, as shown in Fig. 2. 

Summarizing the above considerations, we have that 
for weak, strong, and intermediate system-environment 
coupling strength, maximizing || AH || is expected to max- 
imize the || Aif ||, and as a result, the properties of PS can 
be connected with a maximization of ||Aff ||. Therefore, 
so long as the impact of AH and AK can be quanti- 
tatively captured by their norms ||Ai?|| and ||A2f||, we 
may consider ||Ai?|| only. In particular, by maximiz- 
ing ||AiI||, we expect that the two states \<fi a (t)) and 
\(j)p{t)) in Eqs. ((7]) and ([8]) will evolve most differently, 



thus yielding an approximate method of predicting and 
understanding the emergence of PS. 

As also stressed in the manuscript, since AH is im- 
portant in both equations of ([7]) and ([5]) [—AK contains 
AH in Eq.©], the contribution of the c-number part of 
AH (i.e., H^ a — Hpp) is nontrivial in determining the 
PS. Without this c-number term, maximization of AH 
would incorrectly predict that eigenstates of Hj are the 
PS regardless of the value of e. This interesting feature is 
much different from a pure-dephasing problem. Indeed, 
in a pure-dephasing problem, there would be no popula- 
tion transitions between \cj) a {t)) and \4>p{t)), and hence 
this c-number term in a pure-dephasing problem would 
be of no interest: it only induces a phase difference be- 
tween \<f>p(t)) and \<fip(t)} and therefore it does not cause 
the overlap (<f>p(t)\cf> a (t)) to decay. 

As two consistency checks, we first return to the first 
limiting case of u x 3> s. In this case 

AH ~ H s aa - Hip. _ (34) 

Because eigenstates of H$ correspond to the maximum 
of \H^ a — H§g\, hence, a maximization of ||AiJ||, the 
energy eigenstates of H$ should be the PS. In the second 
limiting case < £, we have 

AH^e{{a\a s z \a)-(P\a s z \P))aj. (35) 

Therefore, the eigenstates of trf and hence the eigenstates 
of the interaction Hamiltonian Hi, can maximize ||Aff ||, 
which is again in agreement with the known result that 
eigenstates of Hi are the PS in this strong coupling case. 
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